Weak ferromagnetism and internal magnetoelectric effect in LiFeP 2 07 
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The magnetic, thermodynamic, and pyroelectric properties of LiFeP2C>7 single crystals are in- 
vestigated with emphasis on the magnetoelectric interaction of the electrical polarization with the 
magnetic order parameter. The magnetic order below Tjv — 27 K is found to be a canted antiferro- 
magnet with a weak ferromagnetic component along the 6— axis. A sharp peak of the pyroelectric 
current at Tjv proves the strong internal magnetoelectric interaction resulting in a sizable polariza- 
tion decrease at the onset of magnetic order. The magnetoelectric effect in external magnetic fields 
combines a linear and a quadratic field dependence below Tjv. Thermal expansion data show a large 
uniaxial magnetoelastic response and prove the existence of strong spin lattice coupling. LiFeP2 07 
is a polar compound with a strong interaction of the magnetic order parameter with the electric 
polarization and the lattice. 

PACS numbers: 75.30.-m, 75.40.Cx, 75.85.+t, 77.70. +a 



I. INTRODUCTION 

The coexistence and mutual interaction of magnetic 
and electric orders in matter has long inspired scien- 
tists because of a wealth of novel physical phenomena 
observed in those systems over the course of more than a 
century. The magnetoelectric effect, describing the inter- 
action and cross correlation of magnetic (electric) orders 
with electric (magnetic) fields, has now been revealed in 
a large number of materials, starting with &2O3 in 1959 
to 1961pi~— The magnetic and electric fields can be ex- 
ternal as well as internal (due to spin or charge orders), 
the coupling can be linear as well as higher order. A re- 
cent review lists a large number of linear magnetoelectric 
compounds that have been discovered since then£ The 
internal magnetoelectric coupling is realized in materials 
with coexisting magnetic and ferroelectric (FE) orders, 
so called multiferroics^— 

Multiferroic materials can be separated into two 
classes. The first class includes compounds in which 
the FE order is improper and the polarization is in- 
duced by an inversion symmetry breaking magnetic or- 
der, for example the transverse spin spiral in TbMnOs, 10 
Ni 3 V 2 8f ii and MnW0 4 ^ The second class includes 
materials that are already polar (FE) before magnetic or- 
der sets in at lower temperature. The classical examples 
are the hexagonal manganites which become ferroelectric 
above room temperature and develop a frustrated antifer- 
romagnetic order below 100 K . 13 ' 14 A significant coupling 
of magnetic and dielectric properties has been revealed 
in those compounds i 15 ' 16 In general, any polar material 
containing magnetic ions with strong interactions lead- 
ing to magnetic long range order at lower temperatures 
is a possible candidate that warrants further studies with 
respect to the the cross coupling of the magnetic order 
parameter and the electrical polarization. 

LiFeP207 crystallizes in a polar structure with the 
monoclinic space group P2i (No 4). The magnetic ion 
Fe 3+ gives rise to antiferromagnetic (AFM) order at 



about 22 (±5) K, as shown in earlier studies of poly- 
crystalline samples and powders i 17 ' 18 The coexistence of 
the magnetic order with the electrical polarization be- 
low the magnetic ordering temperature makes LiFeP207 
an interesting candidate for studying the magnetoelectric 
coupling in electrically polarized materials, similar to the 
multiferrois of the second class, as defined above. The 
possible interaction of the magnetic order with the elec- 
trical polarization in LiFeP20y and the effect of external 
magnetic fields (magnetoelectric effect) are of particular 
interest and will be investigated in this study. 



II. SYNTHESIS AND EXPERIMENTAL 

Large single crystals (centimeter size) of LiFeP2C>7 
have been grown by a top-seeded solution growth method 
using a mixture of polycrystalline LiFeP207, LiH^PCj, 
and NH4H2PO4. Details of the crystal growth are de- 
scribed elsewhere.— The phase purity of the single crys- 
tals was confirmed by powder X-ray diffraction. 

The magnetic properties were measured in a com- 
mercial magnetic property measurement system (MPMS, 
Quantum Design) along different crystallographic orien- 
tations. The effect of the magnetic order and exter- 
nal fields on the electrical polarization was determined 
from integrating the pyroelectric and magnetoelectric 
currents measured upon changing temperature and field, 
respectively. The lattice response to the magnetic phase 
transition was studied via thermal expansion measure- 
ments employing a high-resolution capacitance dilatome- 
ter. The heat capacity was measured using the relaxation 
method implemented in the physical property measure- 
ment system (PPMS, Quantum Design). 

Fig. 1 shows the structure as viewed along the a— axis. 
The almost regular FeC>6 octahedra and the diphosphate 
groups P2O7 form alternating layers stacked along the 
b— axis. The diphosphate groups and the FeOg octahe- 
dra share corners only. The Li ions are located in an off 
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FIG. 1. (Color online) Structure of a LiFeP2C>7 viewed along 
the a-axis. 



center position in tunnels running along the a— axis (Fig. 
1). The only magnetic exchange interaction between two 
neighboring Fe 3+ ions is the super-super exchange via 
two oxygen ions involving the stretched PO4 tetrahedron. 
The relevant magnetic exchange pathways have been dis- 
cussed recently in more detail.— The magnetic order of 
the Fe 3+ spins (S=5/2) was described based on pow- 
der neutron scattering as AFM with the spin orientation 
mainly along the a— axis, although a small component 
along c provided a better fit to the neutron spectra.— 
The magnetic modulation vector, if =(0,0,0), only allows 
for the AFM spin order of the two magnetic ions within 
the unit cell, the magnetic unit cell is therefore identical 
to the crystallographic cell. 



III. RESULTS AND DISCUSSION 

A. Magnetic order below Tjv=27 K 

The single crystals of L1FCP2O7 grown for this study 
permit a more detailed investigation of the magnetic 
properties. In contrast to the earlier reports ; 17 ! 18 the 
magnetic moment M& measured in low fields (10 Oe) 
along the b— axis shows a sharp increase at Tpj ~ 27 
K, indicating a ferromagnetic (FM) moment along this 
direction (Fig. 2a). The data taken after cooling in zero 
magnetic field (ZFC) and cooling in a 10 Oe field (FC) 
used for the measurement of the magnetization are iden- 
tical. This appears unusual since ZFC and FC data of 
a ferromagnet are commonly different because of mag- 
netic domain formation. However, if the magnetically 
ordered state is a single domain state, both data sets 
will be identical and the ferromagnetic M-H hysteresis 
loop should be very sharp. The inset of Fig. 2a shows 
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FIG. 2. (Color online) 6— axis magnetization of LiFeP207. 
(a) M b vs. T at 10 Oe. (b) M 6 vs. H at different tempera- 
tures. Note that the zero-field cooled (ZFC) and field cooled 
(FC) data in (a) perfectly overlap indicating the single do- 
main property of the magnetic order. The inset in (a) shows 
the ac susceptibility near the phase transition. 



the ac-susceptibility, \' b (measured at 117 Hz with a field 
amplitude of 1 Oe), in a narrow temperature range near 
T^r, between 26.7 K and 27.1 K. The sharp peak of x't, 
with a width of only 0.01 K is consistent with the second 
order transition at the critical temperature Tjv=26.86 K. 

The single domain state is further supported by M- 
H measurements shown at different temperatures in Fig. 
2b. The FM hysteresis loops are indeed very sharp, even 
at temperatures as high as 25 K, just below Try. The 
reversal of the FM moment at the coercive field happens 
instantaneously within the data resolution. Some data 
have been collected in field steps as low as 10 Oe and 
no intermediate magnetization values could be detected 
near the coercive field, i.e. the spontaneous FM moment 
flips at once from positive to negative values and vice 
versa, as indicated by the vertical dashed lines in Fig. 2b. 
Once the direction of the FM moment is set, the b— axis 
magnetization further increases linearly with the field up 
to the maximum field of 50 kOe. The magnetic hysteresis 
loops show a small asymmetry the origin of which is not 
clear. It could be related to the complex interaction of 
the FM order parameter M Y with the AFM components 
L x , L z , and the polarization P Y , as discussed in Section 
III D below. 

The magnetic properties measured along the b— axis 
suggest that the true magnetic state of LiFeP 2 07 be- 
low Tjv is a canted antiferromagnet with the spins tilted 
slightly from the a — c plane towards the b— axis, requir- 
ing a revision of the earlier results obtained from pow- 
der sample a 17 ' 18 and theoretical considerations^ This 
canted spin picture is confirmed by comparing the mag- 
netization measured along different crystallographic ori- 
entations. Fig. 3 displays the magnetization data at 
1000 Oe along the a—, b— and c— axes. Mb shows the 
FM response similar to the low-field data, however, M a 
decreases sharply at T^r, as expected for an antiferromag- 
net. M c experiences only a very minute anomaly at T^v 
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and remains nearly constant at lower temperatures. The 
Curie- Weiss extrapolation of the high-temperature data 
reveals a negative Weiss temperature O ~ — 50 K and an 
effective magnetic moment of fJ-eff — 6 hb, close to the 
expected value for spin 5/2. The negative Weiss temper- 
ature proves that the dominant magnetic interactions are 
antiferromagnetic in nature and the FM moment along 
the b— axis most likely originates from the canted AFM 
order. 

It should be noted that the FM component of the 
magnetic order parameter does not violate any sym- 
metry constraints. According to the neutron scatter- 
ing data of LiFeP 2 07 powders, the magnetic modula- 
tion vector is if =(0,0,0) and the magnetic group of if 
can be decomposed into two irreducible representations, 
r=3(r!-|-r2) J£ The basis vectors of the Ti representa- 
tion include two AFM components (L x — — and 
L z = Si — S2), both observed in neutron scattering 
data, and one FM component (M Y = SY + S%). ^\ and 
^2 refer to the spin vectors of the two iron ions in the 
unit cell. M Y allows for the ferromagnetic moment along 
the 6— axis, as observed in our measurements. The fact 
that the FM rise of M fc and the AFM drop of M a happen 
at exactly the same temperature Tjv proves that both 
are components of one and the same magnetic order pa- 
rameter. There remains the question why the b— axis FM 
moment has not been detected in earlier magnetization 17 
and neutron scattering 18 experiments. In both previ- 
ous studies, polycrystalline or powder samples had been 
investigated. Magnetic measurements of randomly ori- 
ented particles can only detect a magnetization averaged 
over all orientations and the FM moment could have es- 
caped the attention. Neutron scattering measurements 
have a finite resolution limit to detect the size of mag- 
netic moments. The spontaneous FM moment measured 
at 5 K is as small as 0.024 ^s/Fe (Fig. 2b). Therefore, 
powder neutron scattering experiments could not detect 
the FM component of the order parameter. 

The single domain property and the smooth increase 
of the FM order parameter at and below T^r strongly 
suggest the second order character of the magnetic phase 
transition. This is further supported by the results of 
heat capacity (C p ) measurements, shown in Fig. 4. 
C P (T) exhibits a very sharp peak at Tjv, similar to the 
A-shaped anomaly of C p that is frequently observed near 
second order phase transitions if critical fluctuations of 
the order parameter are significant. Since the magnetic 
phase transition is continuous, the expansion of the free 
energy with respect to the order parameters will provide 
a mean field description of the magnetic phase transition, 
as discussed in see Section III D. 



B. Pyroelectric and magnetoelectric response of 
LiFeP 2 7 

Because of the polar structure of LiFeP2 07, a macro- 
scopic electrical polarization exists at all temperatures 
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FIG. 3. (Color online) Magnetization M/H of LiFeP 2 7 mea- 
sured at 1000 Oe along the a—, 6— and c— axes. 
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FIG. 4. Heat capacity of LiFeP2 07 near the magnetic phase 
transition. 



below ambient. However, it is not clear whether or not 
the polar distortion disappears at some higher critical 
temperature resulting possibly in a structural and ferro- 
electric transition. Polarization measurements at ambi- 
ent temperature in electric fields up to 24 kV/cm did not 
detect any hysteresis in the P(E) data. 19 High tempera- 
ture data for LiFeP2C>7 are rare. Magnetic data collected 
up to 700 K have not indicated any additional anomaly 
that could be associated with a high-temperature phase 
transition^ Although high-temperature measurements 
of the ionic conductivity of LiFeP2C>7 have shown a 
crossover near 700 K, differential thermal analysis data 
have not shown any notable anomaly that might suggest 
a structural or ferroelectric transition^ 

The single crystals used in this study have been found 
to be good insulators at and below ambient temperatures. 
Therefore, the change of the electrical polarization can be 
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measured using the pyroelectric current method. Current 
data acquired upon increasing and decreasing tempera- 
ture (rate of temperature change: 3 K/min) along the 
axis of polarization (6— axis) are identical in magnitude 
but with opposite sign, proving that the detected current 
is in fact due to the intrinsic polarization change. Fig. 5 
shows the pyroelectric response upon increasing temper- 
ature. The measured current is significant and it proves 
the strong increase of the polarization towards lower tem- 
peratures. The polarization change between ambient and 
low temperatures, AP&(T), is obtained by integrating the 
current and dividing by the contact area and it is shown 
in the inset of Fig. 5. 

The most prominent feature of the pyroelectric current 
I p (T) is its sharp peak at Tat — 27 K, the critical tem- 
perature of the magnetic phase transition, resulting in 
the drop of the polarization below Tjy. This pronounced 
peak implies a very strong coupling of the magnetic or- 
der parameter and the b— axis polarization, i.e. a large 
internal magnetoelectric effect. A closer inspection of 
the pyroelectric data also shows that AP(,(T) reaches its 
maximum (I p =0) at much higher temperatures of about 
2 Tat, which can be understood as the onset of sizable 
magnetic fluctuations and their effect on the lattice. This 
is also consistent with the heat capacity data showing an 
enhancement of C p over the lattice contribution, starting 
at about 50 K (Fig. 4). 

The external magnetoelectric effect (the change of Pf, 
in applied magnetic fields) was studied at different tem- 
peratures and the results are shown in Fig. 6. The b— axis 
polarization increases in longitudinal fields Hf, below the 
magnetic transition temperature as well as above (35 K 
data in Fig. 6a). However, the magnitude of APf,(H) is 
small at low fields but it increases nonlinearly at higher 
magnetic field. It is interesting that APb(H) appears to 
be linear at low fields for temperatures below Tjy. This is 
shown more clearly by plotting the magnetoelectric cur- 
rent I me (measured at 200 Oe/s) in Fig. 6b. I me increases 
linearly with H;, at all temperatures, however, the inter- 
cept with the vertical axis is finite only at temperatures 
below T^r. The dashed lines in Fig. 6b are a linear fit 
to the data and they show the finite intercept for the 5 
K and 15 K data, but not for the 35 K data. Below Tn, 
Ime jumps to a finite value as soon as the field starts to 
increase. This is a clear indication of a linear magne- 
toelectric effect which is superimposed on the quadratic 
magnetoelectric coupling that exists at all temperatures. 

The superposition of linear and quadratic magneto- 
electric effects below Tjv is a consequence of the sponta- 
neous magnetic order which couples linearly to the field 
and to the polarization, as discussed in more detail in 
Section III D below. It should be noted that the po- 
larization change is symmetric with respect to the field 
orientation and data measured at 5 K in negative mag- 
netic fields have been included as stars in Fig. 6a. Simi- 
lar measurements in a transverse magnetic field have not 
shown any significant change of the b— axis polarization. 
The corresponding magnetoelectric current was an order 
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FIG. 5. (Color online) Pyroelectric current and polarization 
change (inset) vs. temperature of LiFeP20y. 
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FIG. 6. (Color online) (a) Magnetoelectric polarization 
change, APt(H), and (b) magnetoelectric current, 7 me , at 
three different temperatures. The bold stars in (a) show the 
polarization change in a negative magnetic field. The dashed 
lines in (b) are a linear fit of the magnetoelectric current data. 

of magnitude smaller and close to the resolution limit of 
the measurement. 



C. Lattice anomalies at the magnetic transition 

The strong effect of the magnetic order below T^ on 
the lattice polarization requires the presence of significant 
spin-lattice coupling resulting in ionic displacements with 
the corresponding change of the electrical polarization, 
as shown in Fig. 5. The macroscopic distortion of the 
lattice below Tjv is revealed through thermal expansion 
measurements. The relative length changes AL(T)/Lq of 
three orthogonal directions, chosen as the b— and c— axes 
of the crystal and the orientation perpendicular to both 
axes, are shown in Fig. 7 at low temperatures. It is 
remarkable that the crystal's dimensions along b— and 
_L (b,c)— directions change very little in passing through 
the magnetic phase transition. However, the crystal ex- 
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FIG. 7. (Color online) Relative length change of b, c, and 
perpendicular (_L) to b and c with reference to 6 K. The inset 
shows the c— axis thermal expansivity «(T). 



hibits a significant contraction of the c— axis with a pro- 
nounced anomaly at Tn- The thermal expansivity of the 
c— axis, a = d(Ac/co)/dT, increases sharply near Tjv 
and displays a A shaped peak, similar to the pyroelectric 
current and the heat capacity. This demonstrates the 
strong interaction of the magnetic order with the lattice. 

It is not clear at this point why the lattice re- 
sponse to the magnetic order is nearly uniaxial along 
the c— direction only (Fig. 7). From the structure of 
LiFeP2C>7 it is obvious that the diphosphate group is par- 
tially oriented along the direction of contraction. It ap- 
pears therefore conceivable to assume that the magnetic 
order mainly distorts the diphosphate group to maximize 
the energy gain of the magnetic system. A decrease of the 
angle P-O-P within the diphosphate group from 129° will 
also result in a contraction of the c— axis. The distortion 
of the PO4 tetrahedra results in a change of the shape of 
the tunnels running along the a— direction. Therefore, a 
displacement of the Li ions in the tunnels can be expected 
resulting in the observed decrease of the polarization be- 
low T^. 

The distortion of the PO4 units below the magnetic 
phase transition will change the super-super exchange 
interactions which depend strongly on the 0-0 distance 



along the exchange path Fe-O-O-Fe^S The observed lat- 
tice contraction (Fig. 7) is therefore a result of a com- 
petition between the gain of magnetic exchange energy 
due to a distortion-induced modification of the super- 
super exchange parameters and the elastic energy of the 
crystal. Theoretical estimates of the exchange coupling 
parameters^ - suggest that the strongest coupling ( J3) is 
found between the two Fe spins ~§ \ and ^2 in the unit 
cell and a sizable coupling (J4) exists between \ and 
S 2 of two neighboring unit cells stacked along the a— axis 
(here we use the same labels for the Ji as in the calcu- 
lation of Refi^ -). The corresponding exchange parame- 
ters are schematically shown in Fig. 8. The contraction 
along the c— axis will reduce the 0-0 distance of these 
two super-super exchange pathways, increase the mag- 
netic exchange constants, and maximize the energy gain 
in the ordered state. While this discussion based on the 
macroscopic expansion anomalies can only be qualitative, 
high-resolution X-ray or neutron scattering experiments 
could resolve the microscopic details of the lattice distor- 
tion and lead to a more fundamental understanding of 
the observed effects. 



D. Landau free energy expansion 

The second order character of the magnetic phase tran- 
sition is well established through the magnetic (Fig. 2a) 
and thermodynamic (Fig. 4) data. Within a mean field 
description, the Landau free energy can therefore be ex- 
panded with respect to the magnetic order parameter. 
The possible coupling of the components of the magnetic 
order parameter and the lattice polarization will be in- 
cluded. All terms in the free energy expansion have to 
be invariant with respect to the symmetry operations of 
the crystal above T^r and the time reversal operation. 
The crystal's space group P2i has only two symmetry el- 
ements, the identity and the screw type operation (180° 
rotation about the b— axis and a translation by b /2). 
Based on previous and current data, the magnetic or- 
der parameter in the little group of if =(0,0,0) is denned 
by the Ti irreducible representation as (L , L z , M Y ). 
The free energy expansion in terms of the magnetic order 
parameter and the electrical polarization is given by: 



I 

F(t, Nt, P*)=F + ^L 2 + -L* + °-M 2 + + a\L x L z + a 2 L x M Y + a 3 L z M Y - aP Y + ^{P Y f 

+ [X 1 L X L Z + X 2 L X M Y + \ 3 L Z M Y + \ 4 (L X ) 2 + \ 5 {L Z ) 2 + X 6 (M Y ) 2 ]P Y + ... (1) 



Equation (1) includes the second and fourth orders of 
the magnetic order parameters it, and all other terms 
up to third order that are invariant under the symme- 
try operations of the crystal and the time reversal. The 
first line represents the magnetic system, limited to the 



components of the order parameter in Ti representation, 
and the energy of the polarized state. Note that the 
linear term in P Y is symmetry allowed in the P2i po- 
lar structure. The a\, cr 2 , and a 3 terms couple different 
components of the magnetic order parameter and they 
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FIG. 8. (Color online) Schematic plot of the most relevant 
super-super exchange interactions, J3 and J4, between differ- 
ent iron moments (only FeOe octahedra are shown).— 

are commonly of relativistic origin. Neutron scattering 
experiments 1 ^ and our magnetization data (Figs. 2, 3) 
show that the major magnetization component in the 
ordered state is L x whereas L z and M Y are compara- 
tively small. This leads us to conclude that the primary 
order parameter is L x and the two other components 
are coupled to L x through the bilinear terms <7\ and 02 , 
resulting in the increase of all three components of the 
magnetic order parameter at the same critical tempera- 
ture, Tjv. This conclusion is further supported by the 
magnetization data of Fig. 3. The a— axis magnetiza- 
tion shows a clear antiferromagnetic response below Tjv 
with a drastic decrease to lower temperatures. However, 
the c— axis response is far more subtle and a significant 
decrease below Tjv is missing, indicating that the L z 
component of the order parameter is most probably in- 
duced by the weak coupling through the \\L X L z term 
in equation (1). The second line in (1) shows the symme- 
try allowed third order terms which couple the magnetic 
order and the polarization. 

Minimizing the free energy (1) with respect to L x , L z , 
M Y , and P Y will define the thermodynamically stable 
state and the magnetic/dielectric orders as functions of 
temperature. However, because of the large number of 
unknown parameters like a, b,..., <7j, and A^, we will not 
attempt to fit a magnetic solution to the experimental 
data. The electrical polarization can be derived as: 

P Y = I [a - \\L X L z - X 2 L X M Y - X 3 L Z M Y 

-\ 4 {L x f-X 5 (L z f-UM Y f] (2) 

Equation (2) describes the dependence of the polariza- 
tion on the magnetic order parameter. The first term, 
a/ P, represents the temperature dependent polarization 
in the paramagnetic state, as obtained by integrating the 
pyroelectric current above Tjy (inset, Fig. 5). The re- 
maining terms in (2) describe the response of the elec- 
trical polarization to the magnetic phase transition re- 
sulting in the peak anomaly of the pyroelectric current 



and the decrease of the polarization Pb below T/v- Note 
that the significant change of the c— axis length below 
Tjv, as shown in Fig. 7, is not considered in equations 
(1) and (2). A more detailed theory that includes the 
magnetoelastic effects on a microscopic level has yet to 
be developed. 

To include the symmetry allowed terms that couple 
the magnetic order parameter to the lattice strain tensor 
£ii, we have to build invariants of the lowest order in 
L , L z , M Y , and . Since the lowest order terms have 
to be bilinear in the components of the magnetic order 
parameter to preserve the time reversal symmetry, they 
include all combinations of L x , L z , and M Y shown in 
the square brackets of equation (1). Al those terms are 
also invariant with respect to the two symmetry elements 
of the space group, the identity and the screw operation. 
Therefore, only space-group invariant components of 
the strain tensor are allowed to couple in first order to 
the bilinear products of the magnetic order parameter 
components. Those elements are e xx , e yy , e zz , and e xz . 
The thermal expansion data shown in Fig. 7 indicate 
that the major lattice distortion is a compressive strain 
along the c-axis, ruling out any significant contributions 
from e xz (shear strain) and e yy (compression of the b- 
axis). The important terms are e xx and e zz which can 
account for the observed c-axis compression below Tjv- 

The magnetoelectric effects shown in Fig. 6 deserve a 
more detailed discussion. The coupling of the magnetic 
order parameter, the magnetic field, and the polarization 
is of interest. Any product of L x , L z , and M Y with the 
external field H Y is invariant under the operations of 
the space group and the time reversal. The polarization 
P Y is also an invariant. Therefore, the magnetoelectric 
effect is described by trilinear terms in the free energy 
expansion, for example E me oc M Y ■ H Y ■ P Y . In the 
high-temperature phase (T>T/v), M Y = \ Y • H Y {x Y — 
AM Y /AH Y is the magnetic susceptibility), and the mag- 
netoelectric response is quadratic in the external field, 



E me oc X Y (H Y ) 2 P Y (3) 

This is indeed observed above Tjv, see for example the 
35 K data in Fig. 6. 

Below Tjy, however, the existence of a spontaneous 
magnetization M s adds an additional term to equation 
(3). With M Y = M s + X Y ■ H Y , 



E me ocM s H Y P Y + X Y (H Y fP Y (4) 

The existence of a spontaneous magnetization below 
Tj\r now gives rise to a linear magnetoelectric response 
in addition to the quadratic effect, in perfect agreement 
with the experimental results presented in Fig. 6. 
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FIG. 9. (Color online) Scaling plot of the magnetization of 
LiFeP2 07. The critical exponent f3 is determined by the slope 
(dashed line). 



E. Critical scaling of the magnetization 

The mean field theory cannot explain the A type 
anomaly of the heat capacity since the critical exponent 
of C p in the Landau mean field theory is zero. Spatial 
fluctuations of the order parameter which become critical 
near Tjv are essential. Evidence for critical behavior is 
derived qualitatively from the A-shaped anomaly of the 
heat capacity and more quantitatively from the tempera- 
ture dependence of the magnetization (Fig. 2a). M(T) is 
expected to follow the relation M (T) = M ■ (1 - T/T N )P 
near the critical temperature. (3=1/2 is the critical ex- 
ponent of the order parameter in a mean field descrip- 
tion; however, a more realistic model description tak- 
ing into account the spatial fluctuations of the order pa- 
rameter leads to different values of j3 depending on the 
model, the spatial dimension, and the specifics of the or- 
der parameter^ The double logarithmic plot in Fig. 9 
shows a significant range where the scaling formula for 
M(T) is fulfilled, as indicated by the dashed line. The 
estimated critical temperature Tjy =26.86 K agrees well 
with the Tjv determined from the sharp ac susceptibility 
peak. The critical exponent /3=0.349 is clearly smaller 
than the mean field value but it is consistent with values 
obtained for typical spin models (Ising or Heisenberg) 
when critical fluctuations are allowed. 22 This shows the 



limits of the mean field theory in describing the second 
order magnetic phase transition. A more advanced de- 
scription has to take into account the spatial fluctuations 
of the order parameter and the critical behavior as well 
as the details of the magnetic order (L x , L z , M Y ) and 
the microscopic interactions. 

IV. SUMMARY 

Single crystals of LiFeP 2 Oy have been investigated 
with respect to the magnetic phase transition at T^r ~ 27 
K and their polarization properties below room temper- 
ature. We show that the magnetic structure below Tjv 
is more complex than previously discussed. The Fe spins 
exhibit a canted antiferromagnetic arrangement with a 
symmetry allowed ferromagnetic component along the 
crystallographic b— axis. Pyroelectric measurements re- 
veal the electrical polarization increasing between room 
temperature and Tjy. A significant internal magneto- 
electric effect is found in form of a sharp A shaped peak 
anomaly of the pyroelectric current at the magnetic tran- 
sition temperature. A similar anomaly of the heat capac- 
ity and the continuous increase of the ferromagnetic com- 
ponent of the magnetic order parameter prove the sec- 
ond order nature of the phase transition. The lattice re- 
sponse to the magnetic transition is uniaxial with a clear 
anomaly of the thermal expansivity measured along the 
c— axis. All experimental data suggest that LiFeP207 is a 
polar crystal with a strong interaction of the magnetic or- 
der parameter and the electrical polarization. The mag- 
netoelectric effect is shown to exist in second order of 
the field above and below the Neel temperature with an 
additional linear component in the magnetically ordered 
state. The trilinear magnetoelectric coupling terms con- 
tributing to the Landau free energy expansion are derived 
from symmetry considerations. 



ACKNOWLEDGMENTS 

We wish to thank Clarina R. dela Cruz for valuable 
discussions. This work is supported in part by the US 
Air Force Office of Scientific Research, the T.L.L. Tem- 
ple Foundation, the J. J. and R. Moores Endowment, 
and the State of Texas through the Texas Center for 
Superconductivity at the University of Houston and at 
LBNL by the US Department of Energy. W. Z. and P. S. 
H. thank the Welch Foundation (Grant E-1457) and the 
Texas Center for Superconductivity for support. 



1 I. E. Dzyaloshinskii, Zh. Exp. Teor. Fiz. 33, 881 (1959) 

2 D. N. Astrov, Zh. Exp. Teor. Fiz. 38, 984 (1960) 

3 D. N. Astrov, Zh. Exp. Teor. Fiz. 40, 1035 (1961) 

4 V. J. Folen, G. T. Rado, and E. W. Stalder, Phys. Rev. 
Lett. 6, 607 (1961) 

5 H. Schmid, in Introduction to complex mediums for op- 
tics and electromagnetics, edited by W. S. Weiglhofer and 



A. Lakhtakia (SPIE Press monograph; PM123, 2003) p. 
167 

6 M. Fiebig, J. Phys. D: Appl. Phys. 38, R123 (2005) 

7 N. A. Spaldin and M. Fiebig, Science 309, 391 (2005) 

8 Y. Tokura, J. Mag. Mag. Mat. 310, 1145 (2007) 

9 S.-W. Cheong and M. Mostovoy, Nature Materials 6, 13 
(2007) 



8 



T. Kimura, T. Goto, H. Shintani, K. Ishizaka, T. Arima, 
and Y. Tokura, Nature (London) 426, 55 (2003) 

11 G. Lawes, A. B. Harris, T. Kimura, N. Rogado, R. J. Cava, 
A. Aharony, O. Entin-Wohlman, T. Yildirim, M. Kenzel- 
mann, C. Broholm, and A. P. Ramirez, Phys. Rev. Lett. 
95, 087205 (2005) 

12 A. H. Arkenbout, T. T. M. Palstra, T. Siegrist, and 
T. Kimura, Phys. Rev. B 74, 184431 (2006) 

13 M. Fiebig, D. Frohlich, K. Kohn, S. Leute, T. Lottermoser, 
V. V. Pavlov, and R. V. Pisarev, Phys. Rev. Lett. 84, 5620 
(2000) 

14 F. Yen, C. dela Cruz, E. Galstyan, Y. Y. Sun, M. Gospodi- 
nov, and C. W. Chu, J. Mater. Research 22, 2163 (2007) 

15 B. Lorenz, A. P. Litvinchuk, M. M. Gospodinov, and C. W. 
Chu, Phys. Rev. Lett. 92, 087204 (2004) 



N. Hur, I. K. Jeong, M. F. Hundley, S. B. Kim, and S.-W. 
Cheong, Phys. Rev. B 79, 134120 (2009) 
D. Riou, N. Nguyen, R. Benloucif, and B. Raveau, Mat. 
Res. Bull. 25, 1363 (1990) 

G. Rousse, J. Rodriguez-Carvajal, C. Wurm, and 
C. Masquelier, Solid State Sciences 4, 973 (2002) 
W. Zhang and P. S. Halasyamani, Cryst. Growth Des. 12, 
2127 (2012) 

M.-H. Whangbo, D. Dai, and H.-J. Koo, Dalton Trans. 
2004, 3019 (2004) 

G. VTthis, Z. Kanepe, A. VTthis, J. Ronis, A. Dindune, and 
A. Lusis, J. Solid State Electrochem. 4, 146 (2000) 
S.-K. Ma, Modern Theory of Critical Phenomena (W. A. 
Benjamin, Inc., 1976) 



